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ABSTRACT 

We present infrared spectroscopy of the Antennae Galaxies (NGC 4038/4039) with NIRSPEC at the 
W. M. Keck Observatory. We imaged the star clusters in the vicinity of the southern nucleus (NGC 4039) 
in 0".39 seeing in K-band using NIRSPEC's slit-viewing camera. The brightest star cluster revealed in 
the near-IR (Mk(0) ~ —17.9) is insignificant optically, b ut coincident with the highest surface brightness 
peak in the mid-IR (12 — 18/im) ISO image presented by Mirabel ct al. (19"98[ ). We obtained high signal- 
to-noise 2.03—2.45 /im spectra of the nucleus and the obscured star cluster at R ~ 1900. 

The cluster is very young (age ~ 4 Myr), massive (M ~ 16 x 10 6 Mq), and compact (density ~ 115 
Mq pc~ 3 within a 32 pc half-light radius), assuming a Salpeter IMF (0.1 — 100 Mq). Its hot stars have 
a radiation field characterized by T c g ~ 39, 000 K, and they ionize a compact H II region with n ~ 10 4 
cm~ 3 . The stars are deeply embedded in gas and dust (Ay ^9 — 10 mag), and their strong FUV field 
powers a clumpy photodissociation region with densities nu ;> 10 5 cm -3 on scales of ~ 200 pc, radiating 
Lh 2 i-o s(i) = 9600 L . 

Subject headings: galaxies: individual (NGC4038/39, Antennae Galaxies) — galaxies: ISM — galaxies: 
starburst — galaxies: star clusters — infrared: galaxies — H II regions 



1. INTRODUCTION 

The Antennae (NGC 4038/4039) are a pair of disk galax- 
ies in an early stage of merging which contain numer- 
ous massive super star clusters (SSCs) alon g their spiral 
arms and around their interaction regi on (Whitmore & 
Schwcizer, 199E; Whitmore et al., 1999] ). The molecular 



~ 1900) K-band spectra of both the obscured cluster and 
the NGC 4039 nucleus. 

2. OBSERVATIONS & DATA REDUCTION 

NIRSPEC is a new facility infrared (0.95 — 5.6/im) spec- 
trometer for the Keck-II t elescope, commission ed during 
April through July, 1999 ( [McLean et al., 1998|) . It has a 
cross-dispersed cryogenic echelle with R ~ 25, 000, and a 
low resolution mode with R ~ 2000. The spectrometer de- 
tector is a 1024 x 1024 InSb ALADDIN focal plane array, 
and the IR slit-viewing camera detector is a 256 x 256 
HgCdTe PICNIC array. 

We observed the Antennae with NIRSPEC during 
the June 1999 commissioning run. Slit-viewing camera 
(SCAM) images at 2 /zm reveal that the mid-IR ISO peak 
is a bright (K = 14.6) compact star cluster located 20". 4 
east and 4". 7 north of the K-band nucleus. This cluster is 
associate d with a faint (V = 23.5) red (V— I = 2.9) source 
(# 80 in Whitmore fe Schwcizer 1995] ) visible with Space 
Telescope (Whitmore & Zhang, private communication). 
We obtained low resolution (R ~ 1900) A ~ 2.03-2.45 ^m 
spectra through a 0".57 x 42" slit at PA=77° located on 
the obscured star cluster and the nucleus of NGC 4039. 
The total integration time on source was 2100 s. 

We dark-subtracted, mean-sky-subtracted, flat-fielded, 
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gas dist ribution peaks at bo th nuclei and in the overlap 
region ( [Stanford ct al., 1990 ), but the gas is not yet un- 
derg oing a global star bu rst typical of more advanced merg- 
ers (Nikola ct al., 1998). Star formation in starbursts ap- 
pears to occur preferentially in SSCs. We chose to observe 
the Antennae because their proximity permits an unusu- 
ally detailed view of the first generation of merger-induced 
SSCs and their influence on the surrounding interstellar 
medium. 

The Infrared Space Observatory (ISO) 12 — 18 /jm im- 
age showed that the hot dust distribution is similar to that 
of the gas, but peaks at an otherwise inconspicuous point 



on th e southern edge of the overlap region (Mirabel et al 



1998). This powerful starburst knot is also a flat-spectrum 



radio continuum source ( Hummel fc van der Hulst, 1986 ) 
and may b e associated with an X-ray source ( Fabbianq 
et al., 1997). We imaged the region around this knot, and 
discovered a bright compact star cluster coincident with 
the mid-IR peak. We obtained moderate-resolution (R 
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and corrected two-dimensional spectra for bad pixels and 
cosmic rays before rectifying the curved order onto a grid 
in which wavelength and slit position are perpendicular. 
We then corrected for residual sky emission and divided by 
a B1.5 standard star spectrum to correct for atmospheric 
absorption. The object spectra were extracted using a 
Gaussian weighting function matched to their strong con- 
tinuua collapsed in wavelength (intrinsic FWHM = 0".84 
for cluster, 0".99 for nucleus) 9 , and then an aperture cor- 
rection was applied to recover the full flux in the continuua. 
Thus we neglected more-extended H2 emission, which has 
maximum FWHM ~ 1".7 in the cluster and ~ 1".2 in the 



nuclei is. We obtained a flux scale by requiring the 2.2 /im 
star flux to equal that corresponding to its K magnitude. 
Reduced spectra are shown in Figures 1 and 2. 
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Figure 1. - NIRSPEC spectrum of the obscured star cluster shows 
nebular and fluorescent H2 emission with a continuum rising toward 
the red. Scaled sky counts are plotted at 0.1 mjy. a>shaped curve 
represents an atmospheric CO2 band at 2.05 fim. 

3. MASSIVE STAR CLUSTER 

The cluster spectrum is characterized by strong emis- 
sion lines 10 and a continuum (detected with SNR ~ 15) 
dominated by the light of hot, blue stars and dust. The 
nebular emission lines are slightly more extended than the 
continuum, and the H2 emission is even more extended. 
This suggests a picture in which hot stars and dust are 
embedded in a giant compact H 11 region surrounded by 
clumpy (see §3.2) clouds of obscuring gas and dust whose 
surfaces are ionized and photodissociated by FUV photons 
escaping from the star cluster. 

For a distance to the An tennae of 19 Mpc (H p=75 
1" 



(Whitmore et al 
-16. 



1998D , we 
We estimate the 



km s" 1 Mpc" 1 , 1"= 93 pc 
find that the cluster has Mr 
screen extinction to the cluster by assuming a range of 
(V— K)p ~ — 1 as e xpected from Starburst99 models 
(|Leitherer et al., 1999]) , and that A K = 0.11 A v QRickc 
L Lcbofsky, 1985|). We find A v = 9 - 10 mag, WEicK 



firmed by our analysis of the H II recombination lines in 
§3.1). We can use the intrinsic brightness along with the 
Lyman continuum flux inferred from the de-reddened Br7 
flux (3.1 x 10~ 14 erg s" 1 cm- 2 ), Q(H+) = 1.0 x 10 53 
photons s _1 , to constrain the cluster mass and age. Us- 
ing instantaneous Starburst99 models we find a total mass 
of ~ 7 x 10 6 M Q (with ~ 2600 O stars) for a Salpeter 
IMF extending from 1 to 100 Mq, and an age of ~ 4 Myr. 
This age is consistent with the lack of photospheric CO 
and metal absorption lines from red supergiants and other 
cool giants, which would begin to contrib ute significantly 
to th e 2 yitm light at an age of ~ 7 Myr (Leitherer et al 



1999Q . The cluster's density is then about 115 M© pc~ J 
for stars of 0.1 — 100 M Q within a half-light radius of 32 
pc. This density is 30 times less than that of the LMC 
SSC, R136 (within a radius of 1.7 pc, assuming a S alpeter 
proportion of low-mass stars) (Hunter et al., 1995). Thus 
the Antennae cluster may be a complex of clusters rather 
than one massive cluster. 
6 
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implies Mk(0) = —17.9, adopting Ak = 1.1 (which is con- 

9 These widths are greater than those measured from the SCAM images 
and rectification errors of order < 1 pixel at the chip edges. 

10 A table of measured line fluxes is available electronically from http://astro.berkeley.edu/~agilbert/antennae 
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Figure 2. - NIRSPEC spectrum of NGC 4039 nucleus shows ex- 
tended collisionally excited H2 emission and a strong stellar contin- 
uum marked by photospheric absorption. No Br7 is present. Scaled 
sky counts are shown at 1.5 mjy. 

3.1. Nebular Emission 

The cluster spectrum features a variety of nebular lines 
that reveal information about the conditions in the ionized 
gas around the cluster, which in turn allows us to constrain 
the effective temperature of the ionizing stars. 

We detected H I Pfund series lines from Pf 19 to Pf 38, 
and display their fluxes relative to that of Br7 in Fig- 
ure 3. The filled symbols give fluxes for the blends Pf 
28+H2 2-1 S(0) and Pf 29+[Fe III]. They fall well above 
the other points, which follow closely the theoretical ex- 
pectation for intensities relative to Br7 (solid curve) with 
no reddening applied, for a gas with n c = 10 4 cm~ 3 and 
T c = 10 4 K flHummer & Storey, 1987| ). Excluding the 
two known blends, the best-fit foreground screen extinc- 
tion is Ak = 1.1 ± 0.3 mag (dashed curve), assuming the 

~ 0."69 and ~ 0".83 (intrinsic), due to the extended line contribution 
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extinction law of Landini et al. (1 



and evaluated at 
2.2 /i m. We consider this an upper limit on Ak because 



a clos i look at the spectrum shows that the points above 
the dashed line in Figure 3 for Pf 22-24 at 2.404, 2.393, 
and 2.383 /im may also be blended or contaminated by 
sky emission, implying a lower Ak and a much better fit 
to the theory. Hence the majority of the extinction to the 
cluster is bypassed by observing it in K band. 
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Figure 3. - Pfund line fluxes relative to Br7 flux (1.05 X 10 — 14 erg 



s cm 
cm -3 , T e 



10 4 K (Hummer & Storey 1987). Filled symbols repre- 
sent lines that are known blends, and the dashed curve shows theo- 
retical fluxes with the best-fit extinction Ak = 1.1 mag. 

The lack of a strong Pfund discontinuity at 2.28 /im in- 
dicates that nebular free-free and bound-free continuum 
is diluted by starlight and dust emission (signaled by the 
rising continuum toward longer A) in the cluster. 

The ratios of [Fe III] lines are nebular density diagnos- 
tics; Table 1 presents observed ratios and theoretical pre- 
dictions of Keenan et al. (1992 ) for emissio n from a colli- 
sionally excited 10 4 K g tabulated by [Luhman et alj 
(1998 ft. The ratios uf [Fe III] 2.140 pm and [Fe HI] 2.243 
/im t o [Fe HI] 2.218 pm are consistent with n c — tQ 3 * 5 — K) 4 - 
cm" 3 . The ratio [Fe III] 2.348 /im/[Fe III] 2.218 /im is 
20% higher than its theoretical val ue, which is roughl y 
constant over all of parameter space ( Keenan et al., 1992j ). 
but [Fe III] 2.348 /im is blended with Pf 29 and subject to 
measurement errors that are larger than the difference in 
extinctions in question (see Figure 3). Even the minimum 
value we infer for this ratio, with Ak — 0, is significantly 
greater than the model prediction. High values of [Fe III] 
2.348 /im/[F e III] 2.218 /im were also found by Luliman| 
et al. (1998) in Orion. This discrepancy may be due to 



blending with another unknown line, or to theoretical er- 
ror; ratios from the latest calculations have an average 
deviation from data of 10% (Keenan et al., 1992). 

He I line ratios can in principal be used to infer nebular 
temperature T e , and are fairly insensitive to n c . However, 
of the three lines we detected, two are not suitable for such 
an analysis: the He I 2.1615+2.1624 /im blend falls on the 
wing of strong Br 7 so its flux has a large (50%) measure- 
ment error, and the strong He I 2.0589 /im line is subject 
to radiative transfer and density effects. 

The He I 2.0589 /im/Br 7 ratio is an indicator of the T e g 
of hot stars in H II regions ( Doyon et al., 1992] ), although 
it is sensitive to nebular conditions such as the relative 



volumes and ionization frac tions of He + a nd H + , geome 
try, de nsity, dustiness, etc. ( [Shields, 1993| ). Dohcrty et al 



(1995j ) studied H and He excitation in a sample of starburst 

For starbursts they found evi- 
(~ 10 2 cm" 2 ) ionized gas from 




galaxies and H II regions, 
dence for high-T e ff, low-n c 
He I 2.0589 /im/Br 7 ratios of 0.22 to 0.64. This is con- 
sistent with giant extended H II regions expected to dom- 
inate the emission-line spectra of typical starbursts. The 
ultra-compact H II regions were characterized by higher 
ratios (0.8—0.9) and higher densities, ~ 10 4 cm -3 . The 
cluster has a flu x ratio of 0.70, a val ue between the two 
object classes of Doherty et al. (1995| ). Assuming the line 
emission is purely nebular, this r atio is consiste nt with a 
high-density (10 4 cm" 3 ) model of [Shields (1993 ), and im- 
plies T e ff ~ 39, 000 K for the assumed model parameters 



This t emperature is similar to that derived by Kunzc et al 
(19960 , - 44,000 K, from mid-IR SWS line observations 
in a large aperture on the overlap region of the Antennae. 
TABLE 1 
Cluster [Fe III] Line Ratios' 1 





Rest 


Observed 


Model Ratio c 


Transition 


A(/im) b 


Ratio 


10 3 


10 4 


10 b 


3 G 3 - 3 H 4 


2.1457 


0.14±0.02 


0.10 


0.17 


0.34 


3 G 5 - 3 H 6 


2.2183 


1.00 


1.00 


1.00 


1.00 


3 G 4 - 3 H 4 


2.2427 


0.28±0.02 


0.26 


0.29 


0.38 


3 G 5 - 3 H 5 


2.3485 


0.80±0.03 d 


0.66 


0.66 


0.66 



Ratios are dereddened fluxes relative to [Fe III] 2.2183 /im, for 



jugar fe Corliss (1985) 



c Models for T C =10 4 K, values of n c in cm -3 ( [Keenan et al., 1992; ). 
d Flux determined by subtracting Pf 29 contribution obtained for 
the best- fit Landini extinction curve with Ak = 1.1 mag. 

The cluster has properties more like those of a compact 
H II region than a diffuse one. It appears to be a young, 
hot, high-density H II region, one of the first to form in 
this part of th e Antennae interaction region (see Habing 



fe Israel 1979 for a review of compact H II regions 
3.2. Molecular Emission 

The spectrum shows evidence for almost pure UV fluo- 
rescence excited by FUV radiation from the O & B stars; 
the strong, vibrationally excited 1—0, 2—1 & 3—2 H2 emis- 
sion has T vib > 6000 K and T rot ~ 970, 1600, and 1800 
K, respectively, and weak higher-v (6—4, 8—6, 9—7) tran- 
sitions are present as well. The H2 lines are extended over 
~ 200 pc, about twice the extent of the continuum and 
nebular line emission, so a significant fraction of the FUV 
(912-1108 A) light escapes from the cluster to heat and 
photodissociate the local molecular ISM. 

We obtained the photodissoc iation region (PDR) mod- 
els of Drainc & Bertoldi (1996) and compared them with 
our data by calculating reduced xl- Models with high 
densities (njj = 10 5 cm -3 ), moderately warm tempera- 
tures (T = 500 to 1500 K at the cloud surface), and high 
FUV fields (Go = 10 3 — 10 5 times the mean interstellar 
field) can reasonably fit the data. Figure 4 shows xl con- 
tours for all models projected onto the njj— Go plane. The 
best-fit Draine & Bertoldi model is n2023b, which has nn 
= 10 5 cm" 3 , T = 900 K, and G = 5000. We fit 22 H 2 
lines, excluding 3—2 S(2) 2.287 /im because it appears to 
be blended with a strong unidentified nebular line at 2.286 
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/ im found in higher- resolution spectra of planetary nebulae 
(Smith ct al., 1981). The weak high-v transitions are all 
under-predicted by this model, and appear to come from 
lower-density gas (nn <J 10 3 — 10 4 cm -3 ) exposed to a 
weaker FUV field (G £ 10 2 - 10 3 ). 

The ortho/para ratio of excited H2 determined from the 
relative column densities in v=l, J=3 and 3=2 inferred 
from 1-0 S(l) and S(0) lines is 1.62±0.07. This is consis- 
tent with the ground state v=0, ,7=1 and J=0 H 2 being 
in LTE with ortho/para ratio of 3 if the FUV absorption 
lines p opulating the non-LTE excit ed states are optically 
thick ( Sternberg k Ncufcld, 1999[ ). Indeed, the best-fit 
PDR models have temperatures that are comparable with 
T rot in the lowest excited states, as well as with the warm 
gas kinetic temperature in the Galactic PDR M16, T = 
930±50 K, measured by Levenson et al. (1999] ). 

If the extent of the H2 emission indicates that the mean- 
free path of a FUV photon is ~ 200 pc, then (nn) = 3 cm -3 
for a Galactic gas-to-dust ratio, while in the PDR(s) nn 
= 10 4 - 10 6 cm . This implies that the molecular gas 
is extremely clumpy, which is consistent with the range of 
densities inferred from the detection of anomalously strong 
v = 8— 6 H2 emission. 

4. NGC 4039 NUCLEUS 

The spectrum of the nucleus of NGC 4039 is marked by 
strong stellar continuum and bright, extended H2 emis- 
sion. Strong photospheric Mg I, Na I, Ca I absorption and 
CO Av = 2 bands indicate that the continuum is domi- 
nated by old giants. The CO band head is stronger than 
that of a M2III, suggesting some contribution from red su- 
pergiants. The absence of Br7 emission implies that star 
formation is currently extinct in the nucleus. Spatially 
extended, collisionally excited H2 emission in the nucleus 
may be excited by SNR shocks from the last generation of 
nuclear star formation, or by merger-induced cloud colli- 
sions. We defer detailed analysis of the nuclear spectrum 
to a later paper. 

5. CONCLUSIONS 

The highest surface brightness mid-IR peak in the ISO 
map of the Antennae Galaxies is a massive (~ 16xl0 6 
M Q ), obscured (A v ~ 9 - 10), young (age ~ 4 Myr) 
star cluster with half-light radius ~ 32 pc, whose strong 
FUV flux excites the surrounding molecular ISM on scales 
of up to 200 pc. The cluster spectrum is dominated by 



extended fluorescently excited H2 emission from clumpy 
PDRs and nebular emission from compact H II regions. 
In contrast, the nearby nucleus of NGC 4039 has a strong 
stellar spectrum dominated by cool stars, where the only 
emission lines are due to shock-excited H2. These obser- 
vations confirm the potential of near-infrared spectroscopy 
for exploration and discovery with the new generation of 
large ground-based telescopes. Our ongoing program of 
NIRSPEC observations promises to reveal a wealth of in- 
formation on the nature of star formation in star clusters. 
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Figure 4. - Comparison of H2 line strengths with PDR models. 
Contours of xt f° r 22 lines projected onto njj— Go plane peak at 
n H ~ 10 5 cm" 3 and G ~ 5000. Model points (+) are for T = 300 
- 2000 K. White + marks best-fit PDR model of Draine & Bertoldi, 
with T = 900 K and xl = 9.3. Contours are 50, 25, 20, 15, 12, 10. 
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